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We investigate the response of a dense monodisperse quasi- two- dimensional (q2D) colloid suspen- 
sion when a particle is dragged by a constant velocity optical trap. Consistent with microrheological 
studies of other geometries, the perturbation induces a leading density wave and trailing wake. We 
also use a hybrid version of Stokesian Dynamics (SD) simulations to parse direct colloid-colloid and 
hydrodynamic interactions. We go on to analyze the underlying individual particle-particle collisions 
in the experimental images. The displacements of particles form chains reminiscent of stress propa- 
gation in sheared granular materials. From these data, we can reconstruct steady- state dipolar- like 
flow patterns that were predicted for dilute suspensions and previously observed in granular analogs 
to our system. The decay of this field differs, however, from point Stokeslet calculations, indicat- 
ing that the non-zero size of the colloids is important. Moreover, there is a pronounced angular 
dependence that corresponds to the surrounding colloid structure, which develops in response to 
the perturbation. Put together, our results show that the response of the complex fluid is highly 
anisotropic owing to the fact that the effects of the perturbation propagate through the structured 
medium via chains of colloid-colloid collisions. 

PACS numbers: 83.50.-v, 82.70.-y, 83.10.-y 



INTRODUCTION 

Provided that one neglects inelastic collisions between 
particles, all fluids should develop hydro dynamic- like dy- 
namics on a sufficiently large time scale. Hence, it is 
reasonable to expect that the same, or at least analo- 
gous, principles will govern all forms of complex fluids. 
To the extent that this is the case, it is important to de- 
termine how the average material response emerges from 
the underlying mechanisms of force transmission and dis- 
sipation. The extreme cases are simple fluids, for which 
the hydrodynamic description is well established [1 , and 
athermal granular bead packs, for which structural char- 
acterizations exist [2]-[5] but descriptions of the dynamics 
remain limited JBH5]- 

Recent theoretical [lOj [11] and experimental [I2j [13] 
advances now allow the study of microscopic force prop- 
agation in complex fluids, such as glasses, colloid suspen- 
sions [10], polymers [14], and granular materials [15]. In 
such experiments p~2j [16] [17] , optical or magnetic tweez- 
ers coupled with optical microscopy are used to micro- 
manipulate a probe particle to deliver a precise and lo- 
calized mechanical perturbation to a system while simul- 
taneously detecting the system response. In this fashion, 
microviscosities in shear thinning and thickening regimes 
have been determined as functions of the speed of the 
probe particle and related to the average structure of the 
fluid around it [H EH CEH] • 

Dense colloid fluids, in which thermal motion and hy- 
drodynamic interactions both play significant roles, lie 
between the liquid and granular regimes. Very recent 
work on sheared colloid fluids [20-22J draws parallels 



to experiments on sheared granular packs [2] while also 
highlighting the complex interplay between granular and 
liquid components. Sheared dense colloid fluids exhibit 
a non-monotonic response with applied shear force, as 
exemplified by shear thinning and shear thickening. The 
recent advances of microrheological methods now make 
it possible to investigate the relationship between these 
non-Newtonian responses and changes in the structure of 
the sheared fluid. We expect changes in structure to play 
an important role in strongly confined fluids because the 
response of a confined fluid to mechanical perturbation 
combines the effects of excluded volume and wall con- 
straints. This reasoning leads to the question of whether 
the viscosity inferred from the force required to push a 
bead inside a microscopic cell via Stokes' Law is appro- 
priate if, for example, the motion of the bead alters the 
local structure of the cell. 

As a step towards establishing the structure and dy- 
namics during the response of a complex fluid to strong 
mechanical perturbation, and examining the role of indi- 
vidual particle-particle collisions, we used localized per- 
turbations to probe the complex properties associated 
with dense colloid fluids. The system that we study is 
a quasi-two dimensional (q2D) monodisperse colloid sus- 
pension. We term it q2D for two reasons. First, by con- 
struction the experimental chamber has a thickness only 
slightly greater than one colloid particle diameter [23] , 
thereby restricting the centers of the colloid particles to 
lie close to a plane but allowing small amplitude devi- 
ations from that plane. Second, the supporting liquid 
in the colloid suspension moves in three dimensions, and 
the functional form, sign, temporal decay, and density de- 



pendence of the effective colloid-colloid interactions are 
affected by the boundary mediated transfer of momen- 
tum and mass in the liquid [24-26 . We note that at the 
densities of our experiments, the system is in the hexatic 
phase. 

First we show that the structure of the fluid around 
one particle driven steadily by an optical trap reprises 
previous observations and calculations [27H29] . Specifi- 
cally, we observe a density wave in front of, and a trailing 
wake following, the probe particle. Surprisingly, our sim- 
ulations suggest that hydrodynamic interactions are not 
required to obtain the density wave and trailing wake, 
although the fluid does influence the size of the wake 
by setting the self-diffusion constant of the colloid parti- 
cles. Second, the angle-resolved pair correlation function 
exhibits a prominent, and heretofore not observed, band- 
ing structure trailing the probe particle. The inter-band 
spacing is consistent with a partial structural rearrange- 
ment of the colloid fluid from hexagonal toward square 
packing. Third, we characterize structurally anisotropic 
displacement chains (excitations of the system) that pre- 
cede the probe particle. We show that the orientation of 
the hexatic order in our dense q2D systems influences the 
formation of these displacement chains. Finally, analysis 
of individual colloid particle displacements reveals a flow 
field around the driven particle with closed circulation 
loops up to 10 particle diameters from the probe parti- 
cle. The high particle density used in our experiment 
induces long range spatial correlations that are evident 
in the observed flow field. Hence, the decay of this field 
deviates from the r -2 behavior predicted for dipolar flow 
using the Stokeslet approximation in q2D systems [24]- 
[26] . Furthermore, the average displacement field shows 
wave-like modulation preceding the probe particle, with 
peaks occurring at integral multiples of the particle di- 
ameter, consistent with the formation of displacement 
chains. Clearly, the fluid is profoundly affected by the 
structure that develops from individual particle-particle 
collisions. 
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FIG. 1. (Color online) Experimental setup, (a) Typical con- 
figuration of q2D monolayer of 2.56 /xm diameter silica beads 
with the probe particle at the center in red with a white 
cross. The field of view is 82/xm x 82/xm and contains ap- 
proximately 1000 particles, (b) Optical trapping setup and 
thickness-controlled sample chamber. 



loop piezo-electric stage between the epi-reflections of a 
reference laser (532 nm) from the top and bottom glass 
surfaces of the sample cell. The chamber walls are coated 
with a non- wetting hydrophobic layer to prevent adsorp- 
tion of colloid particles. Due to the effective repulsion 
between a colloid particle and the walls, the colloid par- 
ticles are confined to the center plane of the chamber and 
form a monolayer with a height distribution of width ±50 
nm. A 810 nm laser beam that is tightly focused with 
a high N.A. objective (Olympus SAPO lOOx oil immer- 
sion) forms an optical trap (Fig. fib) with a force con- 
stant k, = 0.2 pN/nm in the sample plane. This trap 
strength is sufficient to overcome the force threshold (ap- 
proximately 40/c#T/<t) associated with breaking through 
the cage of first neighbors at close packing [10]. All of 
the experiments were carried out in dense colloid suspen- 
sions; the sample was translated at constant velocity past 
the trapped probe particle (i.e., constant dragging). The 
driving speed ranged from 4.9 to 12.6 /xm/s. The total 
field of view (L x L) was 82 x 82 jam 2 and the packing 
fraction, p = Nircr 2 /4L 2 « 0.76. Examination of the 
bond orientational order and the translational order re- 
veals that at this density the system is in the hexatic 
state I301I3T1. 



Experimental System 

The experimental system that we have studied is a 
monodisperse colloid suspension confined between hard 
walls (Fig. []Ji)[30 . The colloid fluid consists of uniform 
size silica spheres (diameter <j= 2.56 ± 0.04 jam; Bangs 
Laboratories) suspended in nanopure water and confined 
to a chamber with 3.2 /am wall separation maintained at 
23° C. The chamber fabrication was adopted from [23] . 
The perpendicular confinement is achieved by modulat- 
ing the hydrostatic pressure inside the sealed chamber 
using a hand pump. The height of the sample cell is de- 
termined by measuring the distance traversed by a closed 



Analysis 

We characterized the strength of the perturbation gen- 
erated by dragging a particle by the Peclet number, 
Pe = <jV/D, where D is the single particle diffusion co- 
efficient at infinite dilution (which we estimate from the 
Stokes-Einstein equation for spheres in 3D), V is the ve- 
locity of the dragged particle and a is the particle diam- 
eter. Physically, Pe is the ratio of amplitudes of driven 
motion to thermal induced Brownian motion for an iso- 
lated particle. Here we reference the Peclet number with 
respect to the diffusion coefficient for an isolated particle 



TABLE I. Summary of experimental parameters 



Drag Speed, V, (pm/s) 


Peclet Number, Pe 


12.6 


190 


8.3 


125 


6.5 


100 


4.9 


74 


Particle Diameter (a) 


2.56 ± 0.04 pm 


Average Packing Fraction (p) 


^0.76 


Field of View (L 2 ) 


82 x 82 pm 2 


Video Frame Rate 


67 Hz 



Ri (i) is a random displacement of particle i satisfying a 
multivariate Gaussian distribution with zero average and 
variance- covariance 



in order to present a consistent measure of the strength 
of the perturbation which is independent of local den- 
sity fluctuations around the probe particle. The particle 
positions and trajectories are extracted from the experi- 
mental data using commercial particle tracking software 
(DiaTrack) . The centers of the particles are determined 
to ~10 nm spatial resolution. Table [I] summarizes the 
important experimental parameters. 



Stokesian Dynamics simulations 

A key characteristic of a colloid suspension is the dis- 
crepancy in size and mass between the colloid particles 
and the molecules of the supporting liquid. The natu- 
ral separation of length and time scales in such a system 
implies that the momentum relaxation of the supporting 
liquid is complete on the time scale of the colloid parti- 
cle motion. As a result, the distribution of particles can 
be described by the Smoluchowski equation. A second 
implication of the natural separation of length and time 
scales is that the coupling between colloid particles can 
be decomposed as (i) a perturbation of the hydrodynamic 
velocity field generated by a single particle, and (ii) the 
corresponding back action of the perturbed velocity field 
on the generating particle. 

Analytical calculation of the effect of the hydrody- 
namic coupling on the colloid diffusive motions in a dense 
system is very complicated [32] [33] . Consequently, we 
use a hybrid version of Stokesian Dynamics (SD) sim- 
ulations to evaluate the contributions of hydrodynamic 
interactions. In such a simulation, the colloid particle 
configuration space trajectories are composed of succes- 
sive displacements, each for a short time step At. The 
influences of particle-particle and particle- wall hydrody- 
namic interactions are incorporated via the values of the 
diffusion tensor D with elements Dij. Specifically, the 
displacement of particle z, Aa^(£), is determined from 
the Langevin equation: 



(Ri(M)Rj(At)) 



2DijAt. 



(2) 



Ax l (At) = J2^ ji &i 



dxj 



J2^At + R t (At). (1) 



Fi is the sum of inter-particle and external forces act- 
ing on particle i and all variables are evaluated at the 
beginning of the time step. 

In principle, for a particular colloid particle configu- 
ration the diffusion coefficient, D^, at time t can be ob- 
tained by solving the Navier-Stokes equation with bound- 
ary conditions supplied by the particle spatial configura- 
tion at the beginning of the time step. We used, in- 
stead, a surrogate procedure that is equivalent. This 
procedure, pioneered by Brady and coworkers [34] [35] , 
captures the effect of the hydrodynamic coupling on the 
diffusive motions of the colloid particles yet greatly sim- 
plifies simulation. In this approach the hydrodynamics 
are captured in the presumed known diagonal and off- 
diagonal diffusion coefficients Da and D^ that appear 
in the Smoluchowski equation. We term our approach 
hybrid Stokesian Dynamics in that the diffusion coeffi- 
cients are derived from experimental data. In an earlier 
study of diffusion in a dense q2D colloid suspension [25] 
we evaluated the center of mass and relative pair diffu- 
sion coefficients, from which the density and separation 
dependences of the diagonal and off-diagonal diffusion 
coefficients can be determined. A fit to those data shows 
that the diagonal and off-diagonal diffusion coefficients 
can be well approximated by Da = a and D^ = f3/r 2 
for particles with center of mass separation r, and that 
both a and /3 are constants depending only on the wall 
separation and the colloid concentration. In our hybrid 
Stokesian Dynamics simulations, in accord with the con- 
ditions appropriate to the experiments reported in this 
paper, these constants are taken to be a = 1 and f3 = 0.3. 

At the beginning of each time step, the diffusion coef- 
ficients Da and D^ were evaluated for the existing par- 
ticle distribution. With these diffusion coefficients, the 
N displacements {Ri} were calculated using a multivari- 
ate normal generator applying the Rotational Method 
described in reference [36 . This method utilizes the fact 
that the k x k variance-covariance matrix I] is positive 
definite and symmetric. Then, if we generate k indepen- 
dent univariate random variables Y = (Yi, Y<2, • • • , Y&), 
we obtain a random sample of X = YP _1 , which 
has a ^-dimensional normal distribution with zero mean 
and variance-covariance matrix X) satisfying P X)P = I, 
where P is the transpose of P. 

It is important to make two further comments concern- 
ing the procedure used to calculate Da and D^. First, 
because they are drawn from experimental data [26j [37] , 
the values of Da and D^, calculated in the fashion de- 
scribed above, incorporate the effects of both far-field 
hydrodynamic and lubrication forces. We note, however, 



that at the densities used in the earlier experiments lu- 
brication forces were found to be insignificant [26l [37] . 
Second, we hypothesize that the particular forms quoted 
for Da and D^, which are valid at packing fraction up 
to 0.49, remain valid at the higher packing fraction of 
the experiments reported in this paper. This idea is con- 
sistent with a result obtained from the point Stokeslet 
representation of the hydrodynamics. Namely, because 
of the no-slip condition at the boundaries that define the 
q2D chamber thickness, the three particle hydrodynamic 
interaction vanishes, although higher order interactions 
do not [25] . We discuss this issue further when comparing 
with the experimental results. 

Our Stokesian Dynamics simulation sample consists 
of N = 85 disks of diameter a in a 10<j x 10<j two- 
dimensional box with periodic boundary conditions. The 
particle-particle interactions were taken to have the con- 
tinuous but nearly hard disk form: 



knT 



C(r* -0.5)- 



(3) 



with r* = r/a, C = 2 x lO" 19 , and 7 = 64. This 
form for the particle-particle interaction has continuous 
derivatives, which simplifies the treatment of forces; it 
has been used previously in simulations of dense colloid 
systems [38]. The system is first brought to equilibrium 
using a Monte Carlo procedure. Then, to simulate the 
experimental arrangement in which a moving laser trap 
is used to exert force on a particle, we drag one of the 
particles in the initially equilibrated sample with a har- 
monic potential whose center moves at a constant ve- 
locity V. The interaction between the dragged parti- 
cle and the center of the moving potential is described 
by Hooke's law, F = — kx, where r is the separation 
of the center of the dragged particle from the center of 
the moving potential and n = 5000/c#T/a 2 . This corre- 
sponds to a trap constant approximately 70 times softer 
than the trap strength used in the experimental system. 
This choice facilitates numerical integration by allowing 
a larger time step but does not significantly impact the 
comparison to experiment because the deviations from 
the trap center are small relative to the particle diame- 
ter in both cases. For the simulations, we use units such 
that a = 1 and Da = 1; time is then measured in terms of 
a 2 /Da. The particle was dragged for 200,000 time steps 
of At = 5 x 10 -5 , after which 100,000 time steps were 
taken to allow the system to relax back to equilibrium. 



RESULTS AND DISCUSSION 
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FIG. 2. (Color online) Pair correlation function g2(r,0) for 
Pe — 125 centered on the probe particle, where r is the inter- 
particle distance and is the angle between the flow direction 
(i.e., dragged particle moving to the left) and the line con- 
necting the particle centers. The grayscale has been chosen 
to highlight the average structure outside the first shell. The 
extreme high (5.5-10) and low (0-2.5) values of the pair cor- 
relation function occur only in the first shell and the wake 
respectively. 




FIG. 3. Pair correlation function #2(7", 0) centered on the 
trapped particle, where r is the inter-particle distance and 
is the angle between the flow direction and the line connecting 
the two particles. The experimental result for Pe = 125 (a) is 
closely reproduced by SD simulations with full hydrodynamic 
interactions (HI), (b). If the HI are turned off and the drag 
speed kept the same (c), the wake diminishes. If the HI are 
turned off and the velocity is increased in order to restore 
aV/Du (d), the amplitude of the wake is comparable to (b). 



Average structure around the probe particle 

We first consider the change in average structure in- 
duced by a probe particle moving at constant velocity 



relative to the bath. Fig. [2] displays the 2D pair correla- 
tion function, g2{r 1 0) 1 defined with respect to the driven 
particle at the origin. The angle defines the direction 



of the projected line between particle centers relative to 
the dragging direction. We use polar binning with ra- 
dial resolution 5r = a/4 and angular resolution 60 = 1°. 
The definition of the pair correlation function in both the 
experiment and the simulation is: 



N 



9(r) 



A{v) i 



53 E *["•-■•*(*)]' 



= 1 2=1 



where r^(t) is the position of particle i at time t and A(r) 
is the area of the bin at position r = (r, 0) relative to the 
probe particle. In other words, for all T frames we bin the 
positions of all N particles using polar coordinates and 
divide the number of counts in each bin by its area. The 
probe, located by definition at the origin, is not counted. 

The trapped particle motion generates a density wave 
that precedes it and a rarefaction in its wake. These fea- 
tures have been previously observed in experiments in 
unconstrained (3D) colloid systems p~3j [27] and in 2D 
simulations [11]. All experimental results shown are for 
Pe = 125, but we observe similar features for 50 < Pe < 
200 (data not shown). In addition to the features al- 
ready discussed (i.e., preceding density wave and trail- 
ing wake), the experimental data shows the emergence of 
clear banding structure trailing the particle. The aver- 
age distance between the peaks of neighboring bands is 
0.88-0.92 a. The absence of comparably prominent and 
sharp banding preceding the particle suggests that the 
motion of the probe induces a partial rearrangement of 
the hexatic fluid. Furthermore, the inter band spacing 
(measured at 3a behind the probe particle) is interme- 
diate between the spacing expected for a square lattice 
(la) and that for a triangular lattice (y/3a/2 « 0.866a). 

To parse the roles of contributing interactions, we per- 
formed simulations for three different constant probe ve- 
locities: crV/Du = 20, 30, and 50. This ratio differs from 
the Peclet number in that Da is not the diffusion co- 
efficient at infinite dilution. Empirically, we found that 
the simulations with aV/Da = 20 gave the best cor- 
respondence to the experiments at Pe = 125, and we 
show results from those simulations unless otherwise in- 
dicated. We focus on this drag speed because the slower 
dragging minimizes the artifacts associated with periodic 
boundary conditions by allowing the system to equili- 
brate within a single periodic unit of the simulation. The 
need to rescale the probe speed likely arises from a num- 
ber of factors: the simulation is truly 2D, the simulation 
box is relatively small in size, and the form of the diffu- 
sion tensor is extrapolated from lower densities. 

Fig. [3] shows the pair correlation function for distance 
ranges well within the periodic boundary conditions. We 
see that the simulations reproduce the leading density 
wave and trailing wake, which suggests that they cap- 
ture the essential physics of the system. To assess the 
importance of the hydrodynamic interactions (HI), we 
compare simulation results with and without HI (Fig. 



|3]3). Without HI we see that the wake still exists, con- 
sistent with the observation of similar patterns in lattice 
Monte Carlo simulations that do not explicitly consider 
frictional forces [39] . Quantitatively, the wake is dramat- 
ically reduced in extent without HI. However, if we in- 
crease the speed of the driven particle to restore aV/Da 
(Fig. [3]i), the pair distribution again resembles closely 
that with full HI. This comparison between simulations 
performed in a consistent fashion clarifies the role of the 
fluid despite potential artifacts arising from the boundary 
conditions. 



Individual displacement chains and hexatic order 

As discussed in the Introduction, our goal is to un- 
derstand how the collective response of the colloid fluid 
arises from the dynamics of individual particle motions. 
In Fig. [4] we plot the displacements of individual colloids 
that move more than a threshold amount. Superimpos- 
ing many such displacements from different times, with 
the probe particle moving over the same range, reveals 
a structure consistent with the density wave and wake 
(black) up to 2 particle diameters from the probe. A 
representative chain of displacements from a small time 
interval of the experiment is highlighted in red. We see 
that it is quite long, spanning 7 particles. Clearly, the 
perturbation propagates via colloid-colloid interactions 
in a highly anisotropic fashion. Displacements chains are 
also observed in our simulations extending up to 4a. Sim- 
ilar behavior has been previously reported in granular 
simulations [40] 

At the densities of our experiments the colloid sus- 
pension is in the hexatic phase. The extent and ori- 
entation of the structural order around a particle can 
be quantified by the 6-fold bond orientation parameter 
^6 — R e {5^7=i ex P(^ ' 6#j)}. Here, 8j is the angle be- 
tween a fixed axis and the line connecting the centers 
of the particle of interest and its j-th nearest neighbor. 
Within regular hexagonal domains, the value of "Q/q is 
constant. Changes in the value of \I>6 correspond to grain 
boundaries and transitions to disordered domains. Fig. 
[5] highlights typical bond orientational configurations ob- 
tained from simulation data under equilibrium and con- 
stant dragging conditions. As expected, in equilibrium 
we observe bond orientational order on length scales com- 
parable to the simulation cell size. For the constant drag- 
ging perturbation, the forced motion of the probe particle 
distorts the hexatic order through repeated cage break- 
ing ultimately leading to the formation of extended grain 
boundaries and disordered regions. 

We would like to determine if the size and orienta- 
tion of the displacement chains is determined by the lo- 
cal structure. Fig. [6] highlights typical situations. The 
sign of ^6 characterizes the orientation of the hexatic 
structure relative to the direction of motion of the probe 




FIG. 4. (Color online) Displacement chains of the bath in- 
duced by a probe particle moving at a constant relative veloc- 
ity (V) in experiment (top) and SD simulations (bottom). A 
collection of particles forms a displacement chain if each par- 
ticle's successive displacements are greater than 0.5d, where d 
is the displacement of the probe particle in the dragged frame 
of reference. The figure depicts a superposition of all displace- 
ment chains observed (black) and a sample of representative 
displacement chains (gray in print or red online). 



particle (compare the orientation of the white hexagons 
with the color scale in Fig. J6|. In each case, the dis- 
placement chain propagates according to the orientation 
of the local structure; displacement chains form inline 
with the probe particle when ^6 > and bifurcate when 
^6 < 0- However, since the size of ordered domains is 
typically less than 10a across, the displacement chains 
can span multiple domains and hence are sensitive to 
grain boundaries. The upper displacement chain in Fig. 
[6^i follows a grain boundary, indicating that both local 
order and defects influence the bath response. The ap- 




FIG. 5. (Color online) Bond orientational order parame- 
ter, ^6, for simulations under equilibrium (left) and constant 
dragging steady state conditions (right). The drag direction 
is to the left and the probe is indicated in red with a white 
cross at (0,0). 




4 -2 

r/CT 



FIG. 6. (Color online) Relation between displacement chains 
and hexatic order. Representative bath particle displacements 
over a time interval of approximately 10s are shown in black 
with the probe particle in red with a white cross at (0,0). 
The color scale indicates the bond orientational order param- 
eter, ^ 6 , calculated for each particle indicated in gray. White 
hexagons illustrate the fact that the sign of ^>6 reports on the 
orientation of the hexatic structure relative to the direction 
of motion of the probe particle. 



proximate angular spread (~ 60°) of the superposition in 
Fig. [4] is consistent with the largest angular spread oc- 
curring when the probe bisects the hexatic ordering, i.e., 
when \I/ 6 = — 6 as shown in Fig. 6a. Orientational effects 
in driving particles through larger ordered domains have 
been investigated previously [4T] . 



Dipolar-like flow 

In the Stokeslet approximation, the perturbation 
caused by the probe particle induces a dipole-like flow 
in the supporting liquid. This prediction is valid for the 
far-field flow in the limit of a dilute suspension of point 
particles. In contrast to these assumptions, in our system 
the particle density is high, the experimental ratio of par- 
ticle diameter to confinement thickness is 0.8, and the far- 



field flow limit is not experimentally accessible owing to 
spatial correlations present in the hexatic phase that are 
comparable to the longest measurable length scales. Fur- 
thermore, given the qualitative connections of the dense 
colloid packing to granular systems [2, 42, 43 j, we expect 
the observed displacement chains (Fig. H| to introduce 
an anisotropic component to the flow field. Nevertheless, 
one might expect that qualitative features of the dipolar 
flow of the liquid to persist. Specifically, we expect the 
closed loop circulation pattern to persist, but not the al- 
gebraic form of the fall off of velocity with the distance 
from the source. We attribute the change in spatial fall 
off of the velocity field to interferences generated by scat- 
tering from the densely packed colloid particles and the 
walls, drawing on our experience with calculations using 
the method of reflections [44]. This method represents 
the influence of the boundaries on the velocity field gen- 
erated by the moving particle as the superposition of the 
velocity components reflected from each boundary en- 
countered by the initially generated velocity field. The 
reflections are considered to be coherent, so the superpo- 
sition of the velocity components has both positive and 
negative interferences. In a system with many reflect- 
ing boundaries, as in the case of densely packed colloids 
studied here, these interferences can greatly change the 
symmetry and spatial falloff of the velocity field. 

We recover the dipolar-like flow shown in Fig. [7] 
through a combination of spatial binning and time av- 
eraging. To construct the field from experimentally mea- 
sured particle trajectories, the frame-to-frame displace- 
ment (A x ,A y ) for each particle was computed in the 
dragged frame of reference: 

(A x , A y ) = (x(t + St) - x{t),y(t + St) - y(t)) , (4) 

where St is the time between frames, and x(t) is the posi- 
tion of particle i at time t. The precise components of the 
drag speed depend on the angular position of the cam- 
era relative to the stage. This introduces an additional 
source of uncertainty, estimated at < 3% in each compo- 
nent, in determining the displacement field components 
in the dragged frame of reference. We choose the drag 
velocity within the uncertainty range that maximizes the 
spatial extent of the dipolar-like flow. This procedure is 
justified by the consistency of the optimization process in 
obtaining far field dipolar-like flow, and the reproducibil- 
ity of the structural form of the angle dependent decay 
in Fig. U 

The displacements were spatially binned based on the 
angle and radial distance relative to the probe particle. 
This was done for every frame of the video, after which 
each bin was time averaged. Only spatial bins that had 
more than 30 samples were used to generate an interpo- 
lated vector field. The threshold of 30 samples was suf- 
ficiently high to ensure that incoherent thermal motions 
were averaged out while ensuring that the whole vector 



field, including regions of vanishing flow, was adequately 
sampled. Fig. [T^i shows the interpolated vector field, de- 
picted by green arrows, superimposed on a line integral 
convolution (LIC) flow visualization [45] which empha- 
sizes the weak but coherent dipolar-like flow observed 
out to ten diameters away from the probe particle. Fig. 
[7)3 shows streamlines (green arrows) and the magnitude 
of the average displacement field (gradient). To the best 
of our knowledge, this is the first experimental observa- 
tion of dipolar-like flow in q2D colloid suspensions, and, 
more generally, the only quantitative analysis of such a 
feature in a complex fluid. 
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FIG. 7. (Color online) Dipolar-like flow fields around a 
probe particle moving at constant velocity for Pe =190. The 
dragging direction is to the left and the field of view has a 
diameter of 20a. The interpolated vector field is obtained 
from the frame to frame displacements of every particle in 
the field of view using a combination of spatial binning and 
time averaging (see text). Fig. [7a shows a sample of the inter- 
polated vectors (green arrows) of the averaged displacement 
field and the associated line integral convolution (gray) indi- 
cating the extent of weak but coherent dipolar flow. Fig. [7J3 
shows streamlines (green arrows) obtained from the interpo- 
lated vector field and the magnitude of the average displace- 
ment field (gradient). 



As already noted, the dipolar-like flow that is gener- 
ated in the supporting liquid by the driven particle is 



expected from a simplified analysis of the hydrodynamic 
behavior of a q2D colloid system wherein the driven par- 
ticle is treated as a force monopole in an incompressible 
liquid confined between two plates [46j [47] . As each col- 
loid particle moves in the plane parallel to the walls it 
changes the momentum and mass density of both the 
colloid particle and the supporting liquid [24]. The mo- 
mentum of the supporting liquid is transferred to the 
confining walls, and is not conserved. The mass of the 
colloid particle and that of the supporting liquid are con- 
served, and the mass perturbation associated with parti- 
cle motion dominates the flow pattern at distances larger 
than the thickness of the sample. To leading order in the 
ratio of particle radius to distance from the particle, the 
mass perturbation can be represented as a mass dipole, 
which in q2D creates a flow velocity proportional to r~ 2 
at large distances from the driven particle. It is this 
functional form that defines the long-range flow pattern 
in q2D confinement. On a time scale short relative to the 
time between colloid-colloid collisions, the colloid parti- 
cles act as reporters of the local fluid velocity, rendering 
the dipolar-like flow observable on a length scale compa- 
rable to multiple colloid-colloid separations. As already 
noted, in a dense colloid suspension we expect deviation 
from the r -2 behavior because of interference with the 
velocity field generated by the zero slip boundary condi- 
tion on each colloid particle with nonzero volume. 

Given the anisotropic structure of both the average 
colloid packing (Fig. [2| and the intermittent colloid 
structure (Fig. [4]), one should expect that colloid dis- 
placements reflect not only the average displacement of 
the underlying fluid, but also the observed structural 
anisotropy. Keeping in mind the roles of (i) colloid-colloid 
interactions, (ii) hydrodynamic screening due to the high 
density, and (iii) nonzero size of the probe particle, we 
determine the decay parameter p in l/r p as the slope 
of a linear fit of log 10 [^6»( r )] as a function of log 10 (r) 
where do(r) is the magnitude of the displacement vector 
of the dipolar-like flow field at a radial distance r from the 
probe particle and fixed angle 6. The fitting region was 
r G (3a, 10a). As shown in Fig. [8^l, the decay parameter 
p exhibits a strong dependence on 0. The observed angle 
dependence may be qualitatively understood as a man- 
ifestation of the nontrivial angle dependent structure of 
the hexatic fluid, as characterized by the pair correlation 
function, and the corresponding colloid-colloid interac- 
tions. We note that some quantitive aspects of Fig.[7]and 
Fig. [8k (e.g., spatial extent of coherent dipolar-like flow 
beyond 5a in Fig. [7] and angular location of minimum 
and quantitative values of the weak decay parameters in 
Fig. [8^l) are sensitive to the components of the drag speed 
within our uncertainty range but the qualitative features 
discussed above are not. 

The modulation of the displacement field preceding 
the probe particle, shown in Fig. [7] exhibits peaks corre- 
sponding to integer particle diameters and a cone consis- 
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FIG. 8. (Color online) Decay parameter of the dipolar-like 
flow field, (a) Angle dependence. The angle corresponding 
to motion of the probe particle is 180 degrees, (b) Typical 
linear fit of logio-logio data for determining the angle de- 
pendent decay parameter. The logio-logio data were fitted to 
f(x) = b — mx and the fitted parameters were m — 1.17 ±0.09 
and b — — 1.40 ± 0.07. The fitted negative slope, m, corre- 
sponds to a single decay parameter quoted in (a). The angle 
shown corresponds to 135 degrees. The uncertainties quoted 
are asymptotic errors in a least squares fit. 



tent with the angular distribution of displacement chains 
shown in Fig. [4] The anisotropic displacement chains, 
which are a direct result of correlated colloid-colloid col- 
lisions, perturb the long range dipolar flow field. Qualita- 
tively similar flow fields were observed in driven granular 
systems [42j [43]. Our simulations, regardless of whether 
HI are included, show dipolar-like flow fields (Fig. [9]) 

7j. We 
analo- 



that are similar to the experimental results (Fig. 
obtain the displacement fields shown in Fig. [9 
gously to the experimental system. In the experiment, 
the displacements (A^AJ were computed using adja- 
cent frames; in the simulation, the displacements were 
computed over the number of simulation time steps re- 
quired for the probe particle to move a distance of a/ 10. 
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FIG. 9. (Color online) Dipolar-like flow fields observed in 
simulations. The field of view has a diameter 10a which spans 
a single simulation cell. The color scale corresponds to the 
magnitude of the displacement field and the streamlines are 
depicted by green arrows. 



In addition to the dipolar-like flow fields, the displace- 
ment field also shows velocity modulation, preceding the 
probe, occurring at integer particle diameters as in the 
experiment. However, the decay is clearly faster than 
in Figs. [7] and [8j Since we obtain similar results even 
for very large systems without HI (data not shown), we 
do not believe this is a boundary artifact. Rather, the 
simulation results suggest that the assumed form for the 
HI, while adequate for elucidating qualitative features of 
the observed dipolar-like flow fields, becomes inappro- 
priate for quantitative analysis of the decay parameters. 
At high densities, an alternative scaling and parameteri- 
zation may be required to capture the interplay between 
granular and liquid components in the colloid fluid, which 
is dominated by near-field effects. 

While a first-principles analysis of experiment and sim- 
ulation is not possible, we can relate the qualitative fea- 
tures observed in Fig. [8] to the phenomenology observed 
in the dipolar-like displacement fields in both experiment 
(Fig. [7]) and simulation (Fig. [9]). As in Fig. |8| we adopt 
the convention that the driving direction is 180° and the 
system is approximately symmetric with respect to the 
drag direction. Consider the following regions: 

• Region I, \0\ E (0°,20°): This region corresponds 
to the wake. Slow decay is observed suggesting that 
the dipolar flow is dominated by the slow structural 
relaxation of the wake. Although the fitting region 
(> 3<j) lies outside the wake (Fig. [2] and Fig. [3]) the 
sampling is poorer in this region than in others. 

• Region II, \0\ G (20°, 60°): This region corresponds 
to dipolar-like circulation currents that flow into 



the wake. Decay parameters with typical values 
p = 1.5 are observed which are between those ex- 
pected for a point force in 3D (p = 1) and the point 
stokeslet limit in q2D (p = 2). 

• Region III, \8\ e (60°, 150°): This region spans the 
location of the vortex. In each of the drag speeds 
consider in the experiment, the minima in angle 
dependent decay at ~ 100° corresponds to the an- 
gular position of the vortex. Comparing the point 
stokeslet approximation and the experimental con- 
ditions, we hypothesize that the non-zero size of 
the probe particle plays a significant role in form- 
ing the vortex. Additionally, we note that both 
simulation and experiment show characteristic an- 
gle dependent anisotropy in the displacement field, 
with the location of the vortex corresponding to a 
steep decline in the displacement magnitude. 

• Region IV, \0\ e (150°, 180°): This region spans 
the velocity modulation of the dipolar-like field di- 
rectly caused by the colloid-colloid collisions. The 
faster decay parameters (p > 2) show that mo- 
mentum transfer through colloid-colloid collisions 
decays more quickly than hydrodynamic interac- 
tions, further suggesting that the long range, co- 
herent dipolar-like flow fields arise from HI. 



CONCLUSIONS 

Our study establishes the phenomenology of response 
to a localized perturbation in a dense q2D colloid sus- 
pension and allows it to be related to individual particle- 
particle collisions. The trapped particle motion generates 
a density wave that precedes it and a void in its wake. 
These ensemble- averaged features have been widely ob- 
served in experiments in unconstrained (3D) colloid sys- 
tems [13l [27] and in 2D simulations [11 . We have demon- 
strated that these features are preserved in a q2D geome- 
try as well and have shown how they arise from individual 
displacement chains. 

Our study provides insight into the range of validity of 
the Smoluchowski equation treatment of the response of 
a colloid suspension to the driven motion of one particle 
PP. This approach neglects multiparticle (more than two) 
hydrodynamic interactions and is thereby limited to low 
density when used to describe a three-dimensional sus- 
pension. It should remain valid to higher densities in our 
case because the three particle hydrodynamic interaction 
vanishes in a q2D colloid suspension at the Stokeslet level 
of approximation, and our simulations did indeed capture 
the leading density wave and trailing wake well. 

That said, the overall deviation of the decay of the 
dipolar-like field from the point Stokeslet approxima- 
tion makes clear the importance of the fact that the 
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colloids are comparable in size to the separation be- 
tween the boundaries. The striking angular dependence 
of the decay reflects the structure of the surrounding 
colloids, which evolves through individual displacement 
chains arising from the perturbation. Incorporating the 
anisotropic propagation of forces and these memory ef- 
fects will be important for theories and simulations that 
seek to describe the response of such complex fluids fully. 
Although doing so will be challenging, success would have 
repercussions not only for colloid suspensions but for a 
much broader range of fluids, including granular ones. 

We will report complementary results from studies of 
relaxation dynamics associated with nonequilibrium ex- 
citations that create density waves, wakes, and persistent 
void areas in a separate publication. 
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